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Abstract. A simple refinement is proposed to the Dibai method for determining black hole masses in 
type-1 thermal AGNs. Comparisons with reverberation mapping black hole masses and host galaxy bulge 
properties suggest that the method is accurate to ±0.15 dex. Contrary to what was thought when the M, 
- a* relationship was first discovered, it does not have a lower dispersion than the M, ~ Li,uige relationship. 
The dispersion in the M, - Lhost relationship for AGNs decreases strongly with increasing M, or Li,uige- 
This is naturally explained as a consequence of the Af , - bulge relationships being the result of averaging due 
to mergers. Simulations show that the decrease in dispersion in the M, or Li,uige relationship with increasing 
mass is in qualitative agreement with being driven by mergers. The large scatter in AGN black hole masses 
at lower masses rules out significant AGN feedback. A non-causal origin of the correlations between black 
holes and bulges explains the frequent lack of supermassive black holes in late-type galaxies, and the lack of 
correlation of black hole mass with pseudo-bulges. 

Keywords: Black hole physics, galaxies; active, quasars; emission lines, galaxies; bulges, galaxies; forma- 
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1 Introduction 



Dibail ([1977|) introduced the estimation of the masses, M,, of supermassiv e black holes (SMBHs 



see 



thermal AGNs (i.e., high accretion rate AGNs seen close to pole-on 
optical luminosity and the FWHM of the H/? line. Because of its ease of use, the 



Antonucci 1993. 20111 



in type-1 
from the 



Dibai method" (also known 

as the "photoionization method" or the "single-epoch method") is by far the most widely used method of S MBH 
mass determination. The assumptions in the Dibai method are discussed in B ochkarev &: Gaskelll (I2009D . For 
a recent review of AGN SMBH determinations see iMarziani fc Sulentid (|201l[ ) . The Dibai method assumes 
that the broad-line region (BLR) is gravitationally bound and that AGNs have similar continuum shapes and 
structures so that the size of the H/3 emitting region can be estimated from the luminosity. For checking these 
assumptions and calibr ating the method, reverberation mapping has been crucial. Cross-correlation of line and 
continuum time series ( Lvutvi fc Cherepashchuk 1972 : Cherepashchuk fc Lvutvi|[l973l : iGaskell fc Sparkd 1986 ) 
yields effective radii of the li ne-emitting r egions, and velocity-resolved reverberation mapping shows that the 
gas is gravitationally bound (|Gaskel]||l988l ). 

A long-standing question in AGN research has been whether the masses of SMBHs are correlated with the 
masses of their host galaxies. Until the mid-1990s it was generally assumed that an SMBH was an incidental 
and perhaps accidental addition to a galaxy - discovering an SMBH in a galax y was considered worthy of a 
press release! A notable exception to this assumption was IZasov fc Dibail (|l970h who discovered a correlation 
between the brightness of AGNs and their host galaxies. iDibai (|l977l) subsequently found an apparent correlation 
between AGN luminosity and M,. Taken together these findings implied that M, and Lhost were correlated. 
The tacit assumption of the 1980s, though, was that any correlation between M, and Lhost was the result of 
selection ef f ects (b ecause only a bright host could be seen in a bright AGN) . This attitude began to change when 
Kormendvl (jl994l) discovered that the masses of inactive SMBHs seemed to be correlated with the luminosity. 
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Lbuige, of the bulge of the host galaxy. Further work f notably iMagorrian et al. 1998 ) verified that there was 
indeed a M, - Lbuige relationshipQ but also showed that there was substantial scatter in it. In trying to 
understand the origin of the relationship it is important to know how much of this scatter is the result of 
observational and modeling errors, and how much is intrinsic. 



2 Improving the Dibai method 



Masses determined from reverberation mapping (jGaskelll 119881 ) have been considered the "gold standard" of 
AGN BH mass determinations because the effective radius is estimated relatively directly. Comparing reverber- 
ation mapping masses, Mrev, with masses from the most widely- used version of the Dibai method {Muitai oc 
FWHM^i"-^) gives a dispersion of ±0.37 dex (see Fig. la). It is k nown that BLR line profile shapes vary from 
object to object and that the shape is a function of line width (see Gaskell 2009bl for a review of the BLR). It is 
therefore reasonable to ask whether the Dibai method can be refined. Gaskelll ( 2010a ) has suggested improving 
the Dib ai method by int roducing an extra line width term. Fig. lb shows Mnew, the masses determined from 
Eq. 1 of iGaskelll (|2010al) . compared with reverberation masses. This gives a relative dispersion of ±0.22 dex. 
An i^-test shows that the improvement is significant at the 99.8% level. 
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Fig. 1. The left frame shows the correlation between masses, Mrev, measured from reverberation mapping and masses, 
Muibai, estimated using from FWHM^L^AGiv- The right frame shows the improvement when Muibai is replaced with 
Mnew, the mass estimated using Eq. 1 of iGaskelli (|2010ai) . H/3 line widths, Lagn, and Mrev are from the compilation of 
IVestergaard fc Petersoiil (|2006l ). 



The ±0.22 dex incl udes errors from both reverbera tion mapping and from using Eq. 1 of IGaskelll (|2010al) . 
From inspection of the Vestergaard fc Peterson ( 2006[ ) compilation of reverberation- mapping black hole esti- 
mates, the rms error in Mrev is ±0.12 dex due to measuring error alone. To this should be added an additional, 
unknown error caused by the off-axis nature of the continuum variability (see Glaskell 2010b). Because of this, 
reverberation mapping gives different time delays for different continuum events. The effect of this on the error 
budge t for Mrev s till ne eds to be evaluated. 

iDennev et al. I (l2009l) have used repeated measurements of the same object to estimate the error in the Dibai 
method due to measuring error. They get an error of ±0.10 dex and, taking into account additional effects, they 
believe the combined error is ±0.12 to 0.16 dex. The rms error in Mnew will be similar. These error estimates 
do not take into account possible random or systematic object-to-object differences in the scaling. 

Since there could be unknown random or systematic sources of error in determination of M, by each method, 
and some of these (such as orientation effects) could be the same for both Alrev and Mnew, it is important to 
have an external check. The correlations between M,, L^uige, and the stellar velocity dispersion, a*, offer such 
a check. 



*The correlation is with bulge luminosity l lKormendv et aLll201lf) , but the Lf^ast measured bv lZasov &: Dibail lll970l ) would have 
been dominated by Lf^^i^^ . 
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3 The M, — Lhost relationship for AGNs 



Bentz et al. ( 20091 ) have presented host galaxy bulge luminosities for most of the AGNs with high-quality 
reverberation mapping mass estimates. Fig. 2a shows the host luminosities compared with masses estimated 
fr om FWHM^L°'^ . Fig. 2b shows the M, - L^uige relationship when the masses are estimated using Eq. 1 
of iGaskelll (|2010al ). As is expected from Fig. 1, the scatter in the M, - Lbuige relationship is reduced. The 
interesting thing is that it is only the scatter at the high-mass end that is reduced. 
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Fig. 2. The improvement in the AGN M, - Lhost relationship when going from using the standard Moibai oc FWHM^L*^'^ 
masses (left-hand frame) to the masses, Mnew, offGaskcU (2010a) (right-hand frame). 



Fig. 3a shows the comparison of residuals from the M, - L^uige and Af, - cr* relationships for the AGNs with 
cr* available. The first thing to note is that the scatter in both axes is comparable. Given that the excitement 
over the discovery of the M, - a ^, relationship was that it was supposed to be significantly tighter than the M, - 
Lbuige relationship (see Fig. 2 of lGebhardt et al.ll2000), this is perhaps surp rising, but making a similar plot for 

bultekin et~aLl (|2009l) (see Fig. 3b) shows that non-BLR 
Lbuiqe and M, - (Th, relationships too. 



the compilation of non-BLR mass determinations in 
galaxies now have comparable scatter about the M, 
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Fig. 3. dlogA/L, the residuals from the M, - Lbuige relationship, versus dlogMcr, the residuals from the M, - cr* 
relationships. The left frame is for the Mnew estimates for the AGNs discussed here and the right frame shows data from 
iGliltekin et all l|2009l ). 

Since M, appears on both axes in Fig. 3, if there are large errors in the mass determinations they will 
introduce a direct correlation into the residuals. No strong correlations can be seen. The slight elongation in 
Fig. 3a corresponds to possible scatter in Mnew of about 0.28 dex. 
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In Fig. 2b it can be seen that the scatter in the M, ~ L^uige relationship declines with increasing mass. Fig. 
4 shows the scatter in the M, - Lhost relationship as a function of Lhost- 
There are several conclusions that can be drawn from Fig. 4. 

1. The scatter f o r the highest luminosity bin (±0.14 dex) is remarkably small. It is comparable to the 
Dennev et al. I (|2009l) estima te of the measuring error. This thus provides good support for the accuracy 



of the approach proposed in Gaskell ( 2010al ). 

2. Since the Mnew determinations seem accurate, the substantial scatter at low luminosities must be intrinsic. 

It is possible that there is some unknown additional source of error causing the increase in dispersion in 
Fig. 4 at low masses, but a similar trend has already been reported for a different sample for the dispersion in 
the M, - cr* relationship (jGaskellll2009afl . Furthermore, the w ide dispersion in M, found from highly-accurate 
maser measurements of yet another sample ( Greene et al. 2010[ ) provides additional strong independent support 
for the large increase in scatter in M,/ Lhost at low masses being real. 
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Fig. 4. The scatter in the AGN M, - Lhost relationship as a function of bulge luminosity. The AGNs have been divided 
equally into four Lhost bins. The dotted line shows the expected n~^^^ reduction in the scatter if a galaxy forms as a 
result of mergers of n units with Gaussian distributions of M, and L stars- The solid line shows the result of simulations 
using a log-normal initial distribution of black hole mass to stellar mass ratio. (Figure adapted from lGaskellll2010ah 



4 The origin of M, - Lhost relationship 



It has lo ng been recognized that mergers will influence correlations between M, and bulge properties. For 
example. iHaehnelt &: KauffmannI ( 20001 ) pointed out that mergers move galaxies along the bulge/ SMBH scaling 
relationships. Pend ( 2007^ ) discussed how mergers might affect the bulge/SMBH scaling relationships and made 
the important point that SMBH-bulge scaling relationships will emerge, independent of the initial conditions 
and detailed physics, simply because of averaging and that "random merging of galaxies that harbor random 
black hole masses tends to strengthen rather than weaken a preexisting, linear, correlation." iGaskeUl (j2010cl ) 
provided the first observational evidence that merging alone is indeed the origin of the Af, - Lhost relationship 
and further proposed that a preexisting correlation is neither necessary nor desirable. If we assume that initially 
M, and the stellar luminosity of each galactic building block have independent normal distributions, then, if 
we merge n building blocks, Lbuige is proportional to n and the scatter in the Al,/ Mh^ige ratio goes down as 
n^^/^ — -^buige- This is shown in Fig. 4 as the dotted line. This is probably too steep a decrease in the scatter^ 
However, a n"-'^/^ decrease is only true for a normal distribution of initial values of Al, and stellar luminosity. 
The actual distributions are almost certainly highly non-Gaussian. Simulations were therefore run starting with 
a more likely log-normal distribution. These gives the solid curve in Fig. 4. 
Several important conclusions follow from these results: 



say "probably" because the right-most point must be somewhat too high because of measuring error, and selection effects (not 
observing very weak AGNs with small SMBHs) could well be lowering the dispersion at low luminosities. 
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1. The M, measurements require a large scatter in M,/Lstars in the initial merging units — a scatter greater 
than that for the least massive galaxies in the sample, where the ±2a scatter is already a factor of ~ 300. 

2. Because of this, detectable SMBHs in So/ Sd galaxies are going to be rare. 

3. Intermediate mass black holes (IMBHs) in globular clusters will also be very rare. 

4. As pointed out in Gaskell ( 2010cl ). the huge scatter in M,/ Lstars in the initial merging units excludes 



fine-tuning of M, or Lhuige or else the subsequent scatter in larger galaxies would be much lower than 
observed. 

Fine-tuning via AGN feedback that significantly reduces scatter at any subsequent stage is also excluded. 
This is not saying, of course, that an AGN has zero effect on its host galaxy — only that it does not 
significantly affect SMBH-bulge relationships. 

The averaging due to mergers will only produce an SMBH-bulge relationship in classical bulges — not 
i n psuedo-bulges which grow through secular processes. This has indeed been found to be the case 
( Kormendv et al.ll201lh . 



The model explains why M, is systematically higher at a given a^, in classical bulges than in pseudo bulges 
(lGrahamll200i Iil]|2008I: iGraham et ai]l201lh , because, in classical bulges, M, is dominated by the highest 
M, in last merger. 



The simulations discussed in iGaskelll (|2010cl) and here are purely for "dry" mergers (those without star 
formation), but SMBHs certainly grow through the accretion of gas as do the stellar populations of galaxies. 
I believe that including star-forming "wet" mergers will not substantially alter the conclusions for two reasons: 
firstly, at the high- mass end we are dealing with galaxies becoming "red and dead". Mergers between them will 
effectively be dry. The second reason why wet mergers will not have a large effect on Fig. 4 is that they are 
most important at the low mass end where the scatter is already large. There are two main possibilities here: 
either the stellar population and the SMBH grow roughly proportionately or they do not. The former gives 
the same result as dry mergers, while the latter can increase the dispersion in the M,/ L,^ta rs ratio. However 



as me ntioned, this dispersion is already very large, so the effect on Fig. 4 will be negligible. iJahnke fc Maccio 
( 201lh have recently performed more elaborate simulations including prescriptio ns for star forma tion, black hole 
growth, and disk-to-bulge conversion, and they have affirmed the conclusions of iGaskelll (|2010a) . 



5 Conclusions 



A slight modification to the Dibai method seems to significantly increase the accuracy of BH mass estimates. 
Using these new mass estimates the scatters about the A/, - Lbuige and M, - cr* relationships are comparable 
(contrary to what was thought when the M, - cr* relationship was discovered). These mass estimates imply that 
there is a strong decrease in the dispersion of the M, - Lhost relationship for higher mass black holes and hosts. 
This is qualitatively consistent with the BH-bulge relationships being solely a consequence of mergers and not 
having an underlying physical cause. The dispersion at the low mass end is too high to allow any significant 
causal relationship between black hole and bulge masses. 
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